Analytical Results for Multifractal Properties of Spectra of 
Quasiperiodic Hamiltonians near the Periodic Chain 



Andreas Riidinger* and Clement Sire 4 " 
Laboratoire de Physique Quantique (URA 505 du CNRS) 
Universite Paul Sabatier, 118 route de Narbonne, 
31062 Toulouse Cedex, France 



Abstract 



The multifractal properties of the electronic spectrum of a general quasiperi- 
odic chain are studied in first order in the quasiperiodic potential strength. 
Analytical expressions for the generalized dimensions are found and are in 
good agreement with numerical simulations. These first order results do not 
depend on the irrational incommensurability. 
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I. INTRODUCTION 



Quasicrystals are characterized by remarkable transport properties (e.g. low conduc- 
tivities increasing with T or disorder) which are generally attributed to the complicated 
interplay between the electrons and the peculiar atomic structure [Q. Some of these have 
been related to the electronic properties found in 2D or 3D theoretical models (critical 
states, spiky singular density of states, anomalous transport...) Q. One dimensional models 
of quasiperiodic structures share all these features, in addition to other properties as the 
occurrence of a zero measure spectrum. Thus, despite their lack of direct physical interest, 
they are worth studying since properties which can only be illustrated numerically in higher 
dimensions, and thus sometimes lack a direct justification, can be calculated very precisely 
in ID. 

A frequently studied one-dimensional quasiperiodic tight-binding Hamiltonian is given 

by 

H = £*i+il*)<* + 1| + Vi\i)(i\ + U\i)(i - 1|, (1) 

i 

where the U 6 {1, p} and Vi G {—A, A} follow quasiperiodic sequences to be specified below. 
The spectra of this kind of Hamiltonians have been proven to be Cantor sets of zero Lebesgue 
measure for all values of the potential 1 — p ^ 0, respectively A / [@|§- Furthermore, 
the density of states (DOS) is purely singular continuous @|§. In order to characterize the 
intricate nature of these spectra, a multifractal analysis of the measure can be elaborated, 
yielding a multifractal spectrum f(a). More precisely, f(a) is the Hausdorff dimension of 
the set of energies in the spectrum with local scaling exponent a, where a characterizes 
the DOS singularities {D(E + AE) — D(E) ~ AE a ). The existence of the generalization 
of such exponents to 2D and 3D can be shown to directly affect the transport properties, 
and a phenomenological transport theory of quasicrystals involve the distribution of scaling 
exponents. Q. In the thermodynamic formalism, introduced by Halsey et ai, generalized 
dimensions D q can be defined 0|| . These also play an important role for the understanding 
of anomalous diffusion properties and the decay of the autocorrelation function encountered 
in quantum dynamics of quasiperiodic systems 0-0. In the case of a quasiperiodic chain 
with the golden mean as incommensurability, a renormalization analysis based on the trace 
map approach leads to the exact minimal and maximal scaling exponents a m i n and a max 
for all values of the quasiperiodic potential, while the fractal dimensions f(a), a m i n < a < 
« max , remain unknown JTD|-|T3J. In the limit of large potential (|1 — p\ 1, |A| 3> 1) the 
renormalization approach of Niu and Nori has been used to obtain f{a) |H>| , |T6f for this 



Fibonacci chain. 

In the present paper we discuss the opposite limit (|1 — p\ -C 1, |A| <C 1), for which a 
perturbation theory approach has been introduced in [|17|] for any quasiperiodic chain. We 
begin by reporting briefly on this approach (II), exploit the results in order to calculate 
the multifractal properties of the energy spectrum (III) and finally compare our theoretical 
relations with numerical data (IV). 

II. PERTURBATION THEORY FOR SMALL POTENTIAL 

We consider the tight binding model, 



2 



(2) 



with the hopping elements tj arranged in a quasiperiodic manner: 

J p for < {ui} < 1 - u, . . 

4 ~ \ 1 for 1 - u < {cot} < 1, 1 j 

where {■} denotes the fractional part and uj 6 (0,1) is an irrational number. The onsite 
potentials are Vi = A if the site z is surrounded by two different bonds, and —A otherwise. In 
the following we will consider rational approximants with p\ hopping elements of strength p 
and qi ones of strength 1 (i.e. ui = 1/(1 + Pi/qi)) and eventually take the limit I — > oo. For 
these approximants with period rii = pi + q\ Bloch theorem can be applied, and one obtains 
rii energy bands. 

For small |A| and |1 — p\ perturbation theory is used, taking as a starting point the 
nonmodulated chain (A = 1 — p = 0) with its eigenvalues, 

Ef(k) = 2cosM^y j = 0,...,n l -l, (4) 

where k represents the Bloch vector. If the quasiperiodic potential is switched on, gaps open 
at the energies Ef ] (0) (j ^ 0, j ^ f (for n { even)) and Ef(ix) (j ^ a±l (for n t odd)). 
For odd values of n\ these energies can be renumbered in a compact form: 

^ = ±2cos(^), J = l,...,^li. (5) 

Since in the limit of large system size, the spectral properties will not depend on whether 
n\ is even or odd, we will restrict ourselves, for sake of simplicity, to odd values of n\. 

In order to facilitate perturbation theory and to index the gaps roughly according to 
their width, Sire and Mosseri [17]] introduced a new numbering, where the atoms of the 
quasiperiodic chain are numbered with respect to their location in perpendicular space rather 
than with respect to their position in physical space. In this basis the originally tridiagonal 
Hamiltonian with quasiperiodic entries becomes a multidiagonal Hamitonian with identical 
entries in each secondary diagonal, except for the two Bloch entries. 

By applying perturbation theory for degenerate states to this Hamitonian it can be shown 
that the band edges read in first order [ |T?1 , |T8[ | : 

Ej ± AEj = 2ecos(7riA)(l - (1 - p)Pi) + A(4fl - 1) 

±^((1 - p) sin(TTjA) - eA sin(27rj'A)), W 

with j G [1, • • ■ , e = ±1 and (3i = pi/ni = 1 — Ui. In principle, this expression correctly 
describes the opening of gaps with index j < 7r 2 /(|A| + |1 — p\). For A = 1 — p = all gaps 
are closed and Ej = 2ecos(7r jpi/n/), which is a mere renumbering of equation (^|). Equation 
(0) is the starting point for our calculation of the multifractal properties. 
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III. MULTIFRACTAL PROPERTIES OF THE ENERGY SPECTRA 



In the thermodynamic formalism the multifractal scaling function r(q) is obtained by 
requiring that the partition function, 



r(r,g) 
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1 1 1 

-Y — 

^ A T 

7 = 1 ^3 



(7) 



be constant for n\ — > oo, where Aj denotes the widths of the individual (nonoverlapping) 
bands. For a periodic chain it is straightforward to obtain that 



9 W(r) 



r + 1 for t < 1, 
2r for r > 1, 



(8) 



so that the multifractal spectrum f(°\a), which is the Legendre transform of r^(q), consists 
of two parts: /^(l) = 1, corresponding to the absolutely continuous component of the 
DOS, and f {0) {l/2) = 0, reflecting the presence of van Hove singularities at the edge of the 
spectrum. 

In the following, our aim is to calculate the corrections to q^°\r), respectively f(°'(a), in 
first order of the quasiperiodic potential. We will separately consider two cases: the onsite 
model (p = 1, A ^ 0) and the off-diagonal model (p ^ 1, A = 0). We begin by treating 
the former model, which is, as we will see, the more interesting one. In order to calculate 
the multifractal properties by means of perturbation theory, we will interpret changes in 
the scaling exponent as logarithmic corrections in the size of the system ni for q(r), in the 



limit of small perturbation (nf 



~ n 



q (0) +Xq (l). 



n(0) 



~ n 



1 + Af^ 1 -* Inn/ + • • ■)). This method 



is reminiscent of what is achieved in the field of critical phenomena in order to calculate 
critical exponents in perturbation theory (e.g. e expansion). 

The bandwidths Aj of the modulated chain can be expressed by 



A? 



A 



3,0 



5A 



3 ' 



where A,- o stands for the corresponding width for the periodic chain: 



A,,o 



cos 



(*m\ _ cos ( <m + 1) 

V ni J \ ni 



2tt 
ni 



sm 



V ni J 



+ 0(- 



n 



2/' 



(9) 



(10) 



and SAj > denotes the decrease of the bandwidth due to the quasiperiodic potential, 
which can be expanded in a formal series, 



6Aj = J2 xl 5A 



i jta(Q 
■j ' 



(11) 



1=1 



For simplicity we take A > 0, so that 5A * > 0; the final result only depends on |A|, which 



can be easily checked. Inserting equation ( |TT| ) m ® and (|?|) we find the partition function 
in first order in A: 



1 + Ar- 



(12) 
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where Vo(t, q) is the partition function of the periodic chain. On the other hand, inserting 
q( t ) = Yl A'gO(r) into the asymptotic relation T(r, q) = nf 9 yields in first order 

r(r,g) = r (r,g)(l + Ag«(r)lnnO. (13) 

Thus the first order correction q^\r) reads 

9 W(r) = r • lim J —= . (14) 

j 

In order to evaluate this expression we remark that the decrease 5Aj = A^o — Aj of each 
bandwidth is the sum of the half-widths of the two adjoining gaps, i.e. in first order: 

X5Af ] = AE h + AE h . (15) 

Since Aj >0 depends smoothly on the location of the band, the sum in the numerator of 
equation (THI) can be rearranged: 



= £ At i- ±|A^|. (16) 



2 



A 



In the diagonal case the half-width of the gap is given by |A£y| = -^j\X sm(2irj so that 
we obtain 



2 1 



^2 -\aJn(icjPi)\ r | cos(7rj$)| 
gW(r) = ^ • lim J- • ^ j . (17) 

- £ |sin(7rj'A)r 

To proceed further we have to distinguish between r < 1 and r > 1. For the first case, 
and provided (3 = lim;_> +00 $ is a generic irrational number (not to well approximated by 
rationals 0]), we can use the Euler-MacLaurin formula to convert the sum in the numerator 
into 



71] — 1 



sin(7rx)| T I cos(Vx) | „,1, 

-dx + C + C(— ), (18) 



Pi x m 



where C is a constant. 

Splitting this integral according to the periodicity of its numerator and replacing the 
sum in the denominator of equation flTTj ) by the corresponding integral, we finally arrive at 
the following expression for q^\r): 



i 

| sin(7rx)| _r | cos(7nr)|da: 

9 {1) (r) = ^- J -^ ri , r<l, (19) 

|sin(7ra)| T dx 



which yields 



^w=^- ( -; + r gL) . r<i. (20, 

Therefore we find the following relation for r(q) in first order in A: 

16 tT(- t+2 ) 

^^' Vr+iirV^r 9 ' 1 - r<L (21) 

Using r(g) = (q — l)D q one obtains for the generalized dimensions 

C « =1 -' A '^ (2- r l?(V) - ?<2 ' ^ 

which remains valid for g — > 2, yielding D2 — 1 — tf|A|. The limit r — ► 1 corresponds to a 
limiting point of the graph f(a), which is given by: 

(23) 
(24) 

The exponent 7 describing the vanishing of the Lebesgue measure or total bandwidth 
(W ~ iM 1 ) satisfies -D_ 7 = and is found to be 7 = 1 — D = ^|A|, in agreement 
with the result given in [pT |. As we have made the assumption that Ag^(r) is small, the 




results for r(q) and D 9 are not valid for large negative values of r or q, where q^'(r) scales 
as \/—t. Therefore it is not possible to obtain a max = D_oo within this approximation. 
To check self-consistency of our results, however, we can set a max equal to the value where 
f{a) = 0. This yields a max = 1 — -^|A 2 , which shows that the decrease of a max < 1 is a 
second order effect in A, in agreement with the results obtained by the trace map for the 



Fibonacci chain [10|-[T3[, and we cannot expect to obtain an exact result for it in this first 
order calculation. We also notice that our result is independent of the irrational number {3. 
In the case r > 1, we find that the limit in equation (|i"7l) vanishes as 1/lnnj for n\ — > 00, 
thus yielding = 0. Therefore the point (a m in = l/2,/(a m in) = 0), which corresponds 

to the van Hove singularities, does not move in first order of the potential. This is also 
consistent with results based on the trace map for the Fibonacci chain. 

We have seen that perturbation theory in first order of the quasiperiodic potential leads 
to multifractal properties for the onsite model (p = 1, A / 0). For the off-diagonal model 
(p 7^ 1, A = 0), however, the half-widths of the gaps \AEj\ = A|l — p\\ sin(7rj/3;)| correlates 
strongly with the corresponding bandwidth of the periodic chain A^q = ^| sm(n j/3i) \. Fol- 
lowing the same line of calculations as for the off-diagonal model, this uniform shrinking of 
the bandwidths leads to 

q^\r) = - 2 .r 1 r < 1. (25) 
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Therefore the corresponding spectrum is a monofractal with all generalized dimensions equal 
to D Q = 1 - £\1 - p\. 

We conclude this paragraph with the remark that it is not the feature diagonal/off- 
diagonal that leads to the different behavior of the two models. Indeed, for the diagonal 
model + Viipi + ipi-i — Eipi with onsite potentials, 



-A for < {ui} <l-u, . . 

A for 1 — uj < {ui} < 1. ^ ' 



the band edges are given by pTTj] 



2A 

Ej ± AEj = 2e cos(ttjA) + A(2/3, - 1) ± — sin(Trj'A). (27) 

The same argument as above shows that in this case the spectrum is monofractal in first 
order in A. 



IV. COMPARISON TO NUMERICAL DATA AND CONCLUSION 



Numerical calculations have been performed to determine the multifractal spectrum. In 
order to speed up the convergence of the numerical data, we calculate the multifractal scaling 
function r(g) by using the more efficient condition, 



r<Q(r,g) 
rP')( r , g ) 



(28) 



rather than by using T^(r, q) = 1 f[9" |. Furthermore, we correct the systematic error due to 
the finite size of the chains by setting 



r«(r,g) 



r('0( T , g ) 



r«(r,g) 



m( T ,q) 



r?(r,g) 



r (n, 

num - 1 - 



7-,g) 



T,g) 



r? } (r,g) 



(29) 



where the subscript "num" refers to the data obtained by numerical diagonalization and 
(r^/To ^)ex is the asymptotically exact value for the undisturbed system: 



rf(r,9) 



Til' 



AO) 



(30) 



This correction procedure is also intended to eliminate the logarithmic contribution to the 
ideal power law ( |30"D expected near (r ~ 1, q ~ 2), which is already present for the periodic 
chain case PD|, and which can strongly affect the numerical determinations of scaling ex- 



ponents pOfl . As we are interested in calculating the first order effects of the quasiperiodic 
potential, we can evaluate r®(r, g)| num according to equation (|T2|), where, of course, SA^ 
has to be replaced by the value A j0 — Aj\ num . 

Finally the corrected value (r^/T^)| C or is inserted into equation ( p8|) for calculating 
r(g), and, afterwards, f(ct) by Legendre transformation. 
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In Fig. [l| the numerical curves f(ct) are compared with the theoretical expression for the 
onsite model and two different values of A (A = 0.001, A = 0.002). The limiting points at 
the left side correspond to r — > 1 and are given by the expressions (^3|) and fl2"4|). The points 
calculated from the numerical data are in good agreement with the theoretical result. The 
difficulty of approaching the limiting point (a c ,f(a c )) is due to the very slow convergence 
of the limit of equation (|T7|) , when r approaches 1. As we have seen, q^\r) shows a 



discontinuity for r = 1, thus preventing good convergence for a finite system even quite far 
from the critical value. This explains why the numerical points end before reaching a = a c . 
However, the numerically accessible a c systematically approaches the theoretical value as 
the system size increases. 

Fig. | shows the dependence of the generalized dimensions D q (q = —5, —2, —1, 0, 1) on 
the quasiperiodic potential A for the onsite model. The lines are the theoretical expressions 
(equation (p2|)). Comparison with the insert, where the corresponding data are plotted for 
the off-diagonal model, shows that the different behavior of the two models that has been 
found theoretically is confirmed by the numerical calculations. 

In conclusion, we have calculated the multifractal spectrum for a general quasiperiodic 
chain in first order in the quasiperiodic potential. We have found analytical expressions for 
the generalized dimensions which are in good agreement with numerical calculations. The 
two models under investigation display qualitatively different behavior, the onsite model 
showing multifractal behavior, the off-diagonal being monofractal. For both models the 
multifractal properties do not depend on the irrational slope in first order in the potential. 
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FIG. 1. The theoretical multifractal spectrum f(a) for the onsite model (p = 1) is plotted for 
A = 0.001 (full line) and A = 0.002 (dashed line) and compared to numerical data for the Fibonacci 
chain with n\ = 987 and ny = 233 (symbols). The theoretical curves are obtained by Legendre 
transformation of equation (21) in the range r G (1,50]. For the numerical data we have chosen 
At = 0.2 for r < 10 and Ar = 5 for r > 10 for sake of clarity. 




strength of quasiperiodic potential X 
FIG. 2. The theoretical dependence of the generalized dimensions D q on the quasiperiodic 

potential A for the onsite model, for q = —5,-2,-1,0,1 (dashed lines from top to bottom) is 

compared to numerical data (symbols). Insert: same as before for the off-diagonal model (A = 0). 
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